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I briefly review the nonrclativistic QCD picture of quarkonium production and its 
■ confrontation with experiment in various production processes. 

■"sj" ' 1 Introduction 

in 
o 

On 



Quarkonium spectroscopy, decay and production has provided us with an in- 
teresting place to test our ideas on QCD ever since charmonium was discovered 
in 1974. Yet, the potential of perturbative QCD (PQCD) to treat production 
and decay has been fully exploited only recentlyEJ in a development comparable 
to that of Heavy Quark Effective Theory for heavy-light mesons. About the 
Qh! same time, experiments measuring quarkonium production at large transverse 

momentum have confronted theorists with surprizingly large cross sections. 
These observations have led_to the understanding that fragmentation □ and 



I _ 

Jv>( , production process. Colorppctet mechanisms were considered in quarkonium 

■ decays already a while agou. They were found to solve the problem of infrared 

divergences in P-wave decays in a systematic way LI. Taking them into account 
also in S-wave production, where they are not required by perturbative con- 
sistency in leading order of a nonrclativistic expansion, opens the promise of a 
quantitative description of quarkonium production. 

In this talk I briefly summarize the concepts underlying the theory of 
quarkonium production. The subsequent survey of the importance of color- 
octet production concentrates on J/ip production. [Of course, most results 
on direct J/ip production generalize to iJj(nS) and T(nS) states.] Because of 
space limitations, the presentation is sometimes sketchy and for many details 
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the reader should consult the original literature, especially to appreciate the 
sources of uncertainty in comparing theory with experiment. 

2 Theory of Quarkonium Production 

2. 1 Factorization 

Inclusive quarkonium production involves two distinct scales. First, a heavy 
quark pair is produced on a distance scale of order 1/mg. Then, the quark 
pair is bound into a quarkonium on a time scale of order of the inverse binding 
energy, r ~ l/(raqu 2 ) (in the quarkonium rest frame), where v is the typical 
velocity of the bound quarks. We assume that v 2 is small (but do not assume 
that the binding force is Coulombic) . The creation process can be computed in 
PQCD and is insensitive to the details of the bound state. The binding process 
can not be computed perturbatively, but long-wavelength gluons responsible 
for binding do not resolve the short-distance production process. The factor- 
ization hypothesis for quarkonium production EJ states that the quarkonium 
production cross section can be written as a sum of short-distance coefficients 
that describe the creation of a QQ pair in a state n multiplied by a process- 
independent matrix element that parameterizes the 'hadronization' of the QQ 
state n into a quarkonium ip plus light hadrons with energies of order vtlqv 2 
in the quarkonium rest frame. Consequently, in a hadron-hadron collision 
A + B — > if) + X, the differential cross section is given by 

d<J = Yl dx 1 dx 2 fi/ A (xi)fj/B(x2) ^2da i+j ^ Q Q [n] (G%), (1) 

i,j n 

where denotes the parton distribution function. The factorization formula 
is diagrammatically represented in Fig. [I] for deeply inelastic scattering. The 
upper part of the diagram represents the matrix element (Of). The hard part 
H is connected to this matrix element by QQ lines plus additional lines, if 
0^ contains more fields. Factorization entails that soft gluons connecting S, 
H and the remnant jet J cancel up to 'higher twist' effects in A, where A 
represents the QCD low-energy scale. Each element in Eq. [l] depends on a 
factorization scale. Since gluon emission changes the QQ state n, a change of 
factorization scale reshuffles contributions between different terms in the sum 
over n and only the sum is physical. The short-distance cross sections da are 
computed by familiar matching: One first calculates the cross section for a 
perturbative QQ state and then subtracts the matrix elements computed in 
this state. 
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Figure 1: Diagrammatic representation of factorization in 7* + A — » ip + X. For the cross 
section the diagram has to be cut. 

A heuristic argument for factorization can be given starting from the in- 
frared fmiteness of open heavy quark production. The above short-distance 
cross sections are obtained by expanding the amplitude squared for open heavy 
quark production in the relative three-momentum of the quarks and by taking 
projections on color and angular momentum states. Since soft gluon emission 
takes one from one state to another, each projection separately is not infrared 
safe. However, by construction the sum in n runs over all states, so that in- 
frared sensitive contributions can always be absorbed in some matrix element. 
In this sense, the sum over all intermediate QQ states restores the inclusiveness 
of open heavy quark production. It is also worth noting that Eq. ^ is valid, 
up to higher twist effects, even if the quarkonium is produced predominantly 
at small transverse momentum with respect to the beam axis, provided one 
integrates over all p t . The transverse momentum distribution is not described 
by Eq. ^ unless p t S> A. A point of concern, however, is, that higher twist 
effects might be of order A/(mQD 2 ), when a soft gluon from a remnant jet 
builds the higher Fock state QQg together with the quark pair. If these terms 
do not cancel, Eq. [l] would be quantitative in hadro-production of quarkonia 
at low p t only for asymptotically large quark masses, since A/(mQV 2 ) ~ 1 both 
for charmonium and bottomonium. 

2.2 NRQCD and Velocity Scaling 

The matrix elements (0%) contain all interactions of the nonrelativistic QQ 
pair (in its rest frame) with degrees of freedom with low momentum compared 
to m,Q. These interactions are accurately described by an effective field the- 



3 



ory called nDnrelativistic QCD (NRQCD). The matrix elements are denned in 
NRQCD asH 

(°n > = EE <°l *W \H(X)X){H(X)X\ ^n' nX |0>, (2) 

X A 

where the sum is over all polarizations A and light hadrons X, and ip, x are 
two-spinor fields. [Matrix elements with additional gluon fields are suppressed 
in v 2 and will not be considered.] Typically, the 'kernels' K n , n' n specify the 
color, spin and orbital angular momentum state of the quark-antiquark pair. 

Without an additional organizing principle, there would be too many ma- 
trix elements to make the theory predictive. The NRQCD Lagrangian informs 
us about the coupling of soft gluons to the QQ pair. In particular, spin sym- 
metry holds to leading order in v 2 and reduces the number of independent 
matrix elements considerably. There is no flavor symmetry as in Heavy Quark 
Effective Theory, because the kinetic energy term is part of the leading order 
Lagrangian. Quarkonium spectroscopy tells us that the kinetic energy is ap- 
proximately constant in the range of reduced masses that comprises cc and bb 
states. Thus, v 2 ~ I/tuq in the range of interest, while for very large quark 
masses v 2 ~ 1/ In 2 toq. Furthermore, the overlap of the final state HX with a 
QQ state n can be estimated from a multipole expansion, which allows us to 
drop some operators that acquire an additional suppression compared to the 
scaling in v 2 of the kernels themselves. The resulting 'velocity scaling rules' 
are summarized in Refs.ulJ. The double expansion of a quarkonium production 
cross section in a s and v 2 (neglecting higher twist corrections) is now complete. 

As a result two matrix elements, (Of c0 ( 3 Po)) and (0g cO ( 3 Si)) are needed 
to describe the production of all three P-wave states at leading order in v 2 . 
The notation refers to the kernels in Eq. The subscript denotes the color 
state and the angular momentum state is written in spectroscopic notation. 
The first matrix element reduces, to leading order in v 2 , to the familiar deriva- 
tive of the wavefunction at the origin. The second color-octet term absorbs 
the IR seastive regions that would otherwise appear in the short-distance cross 
section tffl. At leading order in v 2 , J ftp production is described by the single 
parameter (0(^( 3 Si)}. Because of charge conjugation, the gluon-gluon fusion 
short-distance cross section which multiplies this matrix element is suppressed 
by a s and is proportional to a 3 . Since v 2 ~ 0.25 — 0.3 is not very small for 
charmonium, higher order corrections in v 2 can be important, if they arise at 
lower order in a s . According to the velocity scaling rules three color-octet ma- 
trix elements - {0^( 3 Si)), (O^^So)), {O^( 3 P )) ~ contribute at order 
a 2 v 4 (powers of v 2 are counted relative to the leading order contribution), so 
that J/tp production is described by four nonperturbative parameters. Con- 
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tributions of order a^v 2 also exist and can be important in specific regions of 
phase space. Because Xd production at low tranverse momentum is also sup- 
pressed by a s compared to Xco and Xc2, higher order corrections in v 2 would 
also be important for Xci production. 



2.3 Fragmentation 

We now consider the transverse momentum distribution da /dp 2 . The leading 
order contributions (in a s and v 2 ) decrease as (P-waves) or 1/pf (S- 

waves), when pt is large compared to 2ttiq. This steep decrease is a penalty 
for preforming the quarkonium state at very small distances 1/pt rather than 
1/mq and is not what would be expected from a high-energy cross section in 
QCD. For A + B — > H + X, where H can be any hadron, we expect scal- 
ing: When the cms energy and p t are large compared to all hadron masses, 
the short-distance cross section da /dp 2 scales as 1/pf on dimensional grounds, 
modulo logarithms of p t and higher twist corrections of order mu/pt and A/p t . 
Moreover, the leading twist cross section can be written as a convolution of 
distribution functions, a short-distance cross section and a fragmentation func- 
tion. From their p t -behaviour we deduce that the leading order (in a 3 and v 2 ) 
J/tp production mechanisms arc higher twist at large pt (but calculable, since 
mn is large compared to A). In general, fragmentation functions remain uncal- 
culable. If H is a quarkonium, however, the dependence on the energy fraction 
z can be calculatedoB, because the quark mass tuq 3> A provides another large 
mass scale. Thus, a parton fragments first into a QQ pair, which subsequently 
hadronizes. Since the hadronization of the QQ pair takes place by emission of 
gluons with momenta of order mQV 2 in the quarkonium rest frame, the energy 
fraction of the quarkonium relative to the fragmenting parton, differs from that 
of the QQ pair only by an amount Sz ~ v 2 <C 1. As a result, the fragmentation 
functions are expressed as a sum over perturbatively calculable, z-dependent 
coefficient functions that describe the fragmentation process i — > QQ[rt] multi- 
plied by the same z-independent 'hadronization matrix elements' encountered 
earlier. 

At large pt quarkonium production depends on three small parameters: 
a s /7r ~ 0.1, v 2 r~j 0.25 — 0.3, Aitiq/p 2 (numbers for charmonium). The leading 
contributions for octet vs. singlet and fragmentation vs. non-fragmentation 
production of J/tp are shown in Tab. [l]. The scaling with v 2 and Ajn 2 Jp\ is 
measured relative to the color-singlet non-fragmentation term. We see that at 
Pt ~ 10 GeV gluon fragmentation into a color-octet quark pair dominatesu all 
other mechanisms by at least a factor ten. 
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Table 1: Parametric dependence of various J/ip production mechanisms in hadron-hadron 

collisions. 





color singlet 


color octet 


Pt ~o 




2 4 


Pt > 2m Q 


a 5 s (p?/(4m^)) 2 


a^ 4 (^/(4m^)) 2 



2.4 Quarkonium Polarization 

Qnar koniniTi p olarization in the NRQCD formalism has been considered in 
Refs. lEJO'Ej. A novel feature as compared to unobserved polarization is 
that the short-distance cross sections can not be determined from a match- 
ing calculation that involves only amplitudes squared of QQ states with defi- 
nite angular momentum. The final quarkonium state can be reached through 
quark-antiquark pairs in various spin and orbital angular momentum states, 
which are coherently produced, so that interference between different interme- 
diate states occursBliilO. Because of parity and charge conjugation symmetry, 
intermediate states with different spin S and angular momentum L can not 
interfere, so that for J/ip production interference occurs for intermediate Pj- 
states only. The interference terms are crucial to obtain a.factorized expression 
for fragmentation functions into polarized ip(nS) states El. 

The amplitude projections that determine the short-distance cross section 
can be found as follows. The cross section can be written as (see Fig. |l|) 

^-H ai;bj .S%, (3) 

where the indices ij and ab refer to spin and orbital angular momentum in a 
Cartesian basis L a Si (a,i = 1,2,3 = x,y,z). For QQ with spin and orbital 
angular momentum one, the soft part is given by 

S%=J2<°\x^iT A (- l -D^m(X)X){J/^ 

(4) 

where helicity A is fixed. To proceed one writes down the most general ten- 
sor decomposition compatible with rotational invariance, parity and charge 
conjugation. In general, this introduces a substantial number of new nonper- 
turbative parameters. To evaluate the matrix element at leading order in v 2 , 
we may use spin symmetry. Spin symmetry implies that the spin of the J/ip 
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is aligned with the spin of the cc pair, so S ai : b - oc e l (A)e J (A). One then finds 
S^ bj = {0 J 8 /,p ( 3 Po))6 ab e i *(X)e^X). (5) 

All other possible Lorentz structures are suppressed by v 2 . Moreover, the 
single surviving structure is related to a matrix eLsment that appeared already 
in unpolarized production. This result is generalBr 1 !! 12 !: At order a^ 4 , no new 
matrix elements are required to describe J /if) polarization. 

The decomposition Eq. || tells us that to calculate the polarized pro- 
duction rate we should project the hard scattering amplitude onto states 
with definite S z = X and L z , square the amplitude, and then sum over L z 
Q2l e a{L z )eb{L z ) = S a b in the rest frame). The soft part is diagonal in the 
L Z S Z basis. It is straightforward to transform to the more conventional JJ Z 
basis. Since J z = L z + S z , there is no interference between intermediate states 
with different J z . We write, with obvious notation, 

*n~ £ //,, :/■./■ -.s-,^ (6) 

JJ Z ;J'J' Z 

and obtain, using Eq. [|, 

S J X Lj-J L = (O^CPo)) £<1M; l\\JJ z ){J'J' z \\M; 1A) , (7) 

M 

which is diagonal in (J J Z ){J' J' z ) only after summation over A (unpolarized 
production). In general, the off-diagonal matrix elements cause interference 
of the following JJ Z states: 00 with 20, 11 with 21 and 1(-1) with 2(-l). 
This particular pattern of interference is a consequence of spin symmetry. In 
general, all states interfere, and diagonality in the L Z S Z basis is also lost. 

3 Confrontation with Experiment 

The interest in color-octet production mechanisms was ignited by the possibil- 
ity to explain the large if)' production cross section at the Tevatron by gluon 
fragmentation into a color-octet cc pair. Over the past year most J/tp produc- 
tion processes have been reanalyzed with color-octet mechanisms taken into 
account. In this section I intend to give a short summary. References to ex- 
perimental results can be found in the quoted papers. The letter denotes 
J/tjj and if)' collectively. 
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3. 1 Large pt 



Quarkonium production at large transverse momentum with respect to the 
beam axis can be measured at colliders, most recently in pp collisions at the 
Tevatron, where the indirect contribution from B decays can be removed. The 
accessible range of pt is 5 GeV < pt < 20 GeV. In this range fragmentation 
dominates. The J /ip production data, which includes feed-down from P-wave 
states could be described E3 - within theoretical and experimental errors - by 
fragmentation into P-wave states together with color-singlet gluon fragmenta- 
tion into J/ip- The ip' cross section, however, is under-predicted by a factor 
30 by color singlet fragmentation. The same deficit was found for direct J/tp 
production, after the contributions from higher charmonium states could be 
separated. The dominant S-wave production mechanism was still missed. 

The discrepancy can be explained i0 by taking into account that a gluon 
can fragment into a color octet cc pair in a 3 Si state, which decays into a ip state 
by emission of two gluons with momenta of order m c v 2 in the quarkonium rest 
frame. The suppression factor w 4 is more than compensated by the enhanced 
short-distance coefficient, see Tab. [j] and Sect. 2.3. At lowest order in a s 
only (Og ( 3 Si)) is probed by gluon fragmentation and can be determined by 
fitting the Tevatron data. The resulting (Of ( 3 Si)) / (Of ( 3 Si)) - 1/100 is 
actually smaller than v 4 ~ 1/10 as expected from velocity scaling. Apart from 
numerical factors in the short-distance coefficients, this could partly be due to 
the fact that the emission of soft gluons is kinematically not accounted for in the 
octet fragmentation function. Soft gluon emission softens the fragmentation 
function by smearing the delta-function over a region Sz ~ v 2 in longitudinal 
momentum fraction. Since the short-distance cross section falls off like 1/pf, 
one roughly probes the fourth moment of the fragmentation function. A softer 
fragmentation function would then require a larger matrix element to fit the 
data. 

The theoretical predicition was extended to moderate pt ~ 2m c (or slightly 
larger) in Rcf. £3, where contributions suppressed by Am 2 /p 2 at large pt have 
been kept. In this intermediate region, all three color-octet matrix elements 
relevant for ip production are equally significant, but should be suppressed by 
v A compared to the color-singlet cross section, since no compensating factors 
of 7r/a s are at work. Surprizingly, due to numerical enhancements of the am- 
plitudes, the color-octet contributions still dominate and seem to be required 
to describe the data. The best fit yields a value of l£>g^( 3 Si)) which is more 
than a factor two smaller than the one extracted L3 from large p t , while the 
fitted combination of (O^^Pq)) and (O^^Sq)) is rather large, especially 
in comparison with the fits for t/j'. Whether this difference between J/tp and 
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ip' is an artefact or physical effect remains to be decided. 

As noted in Ref.Ej, the above scenario implies, that the directly produced 
ip mesons are almost completely transversely polarized at large pt- The frag- 
menting gluon is transversely polarized and transfers its polarization to the cc 
pair in the 3 Si state. Because of spin symmetry, the transverse polarization 
stays intact in the emission of two soft gluons. Polarization would be mea- 
sured in the angular distribution of leptonic tp decay, dT/dcos 9 oc 1 + a cos 2 6, 
where 9 denotes the angle between the lepton three-momentum in the tp res t 
frame and the ip three-momentum in the lab frame. Higher order corrections 
to the fragmentation function., can result in longitudinally polarized ip, but 
were found to be rather small!. For (O^( 3 P ))/(m 2 (O^( 3 5i))) < 2, one still 
has a > 0.64. Spin-symmetry breaking corrections are suppressed by v and 
introduce an uncertainty of 0.1 in a. These corrections to a = 1 persist at 
large pt- For moderate pt corrections that vanish as 4m 2 . jp\ exist, but have 
not yet been calculated. Since other corrections are small, these could well 
dominate over most of the pt-range accessible at the Tevatron. So far, the 
angular distribution remains unmeasured. There is a price in statistics to pay, 
but understanding the angular dependence of the detector acceptance is also 
non-trivial. Ideally, a measurement of the pt-dependence of a will provide a 
decisive self-consistency test of the octet production picture. 

At HERA large-pt production of J/tp is probed in jp collisions. In this 
case gluon fragmentation plays no role, since at leading order in a s only "/q 
collisions contribute to this mechanism. Color singlet charm fragmentation, 
induced by photon-gluon fusion, dominates EZI at large pt- 



3.2 Small p t : Fixed Target Experiments 

Fixed target experiments have for a long time been the most profuse source 
of quarkonium production data. Since no cut on p t is usually imposed, the 
production cross section is dominated by quarkonia with transverse momenta 
of about 1 GeV. The cms energies range up to y/s = 40 GcV. The most recent 
comparisons of color-singlet production mechanisms_with data are documented 
in Refs.EjO. Summarizing the conclusions of Ref.E3, the color singlet mecha- 
nisms (a) do not account for the over-all normalization of the total cross section 
very well, (b) yield too low a fraction of directly produced J/tp, (c) transverse 
rather than no polarization of J/tp and ip and (d) predict far too few \ci 
states in comparison to x«s h these failures in mind, color-octet mecha- 

nisms have been analyzed liircJ'EJ. [For reasons explained in Ref.liil, there are 
substantial differences in numerics and conclusions between these three anal- 
yses. The subsequent presentation adheres to Ref. OJ.] The results are not 
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entirely encouraging, but not uninteresting cither. 

Because of charge conjugation, color-singlet -0 production is suppressed 
by a s . Consequently, octet mechanisms scale as (Tr/a s ) v 4 ~ 1, see Tab. [jj 
Since the color-singlet amplitude vanishes at cc threshold while the octet 
amplitude does not, and since the cross section is enhanced at threshold 
by the a>behaviour of the gluon distribution, color-octet production actu- 
ally dominates the cross section. At leading order, the combination Ag(ip) = 
(Ot^So)) +7/m 2 c (Ot{ 3 P )) is probed. At this order (0^( 3 5i)) arises only in 
qq annihilation, which is numerically small compared to gluon-gluon fusion in 
the energy range considered. 

Starting with total unpolarized cross sections, one has three observables 
- crj/ip, cr^< and the direct J/tp cross section CjM' wn i cn does not include 
feed-down from higher charmonium states - and two octet hadronization pa- 
rameters - Ag(J/i/)), Ag(V>') - to fit. [If these were accurately known from other 
processes, no free parameter would remain.] Such a fit is possible with matrix 
elements consistent with their expected size from velocity scaling. The energy 
dependence of the total cross section is in agreement with the data, although 
this can not be considered very significant. A fit to <Jj/^ and cr^i then pre- 
dicts a direct J/tfj production fraction of about 60%, in good agreement with 
experiment. With respect to the failures (a) and (b) above, this could be con- 
sidered as strong evidence for important color-octet contributions in direct ip 
production. The uncertainties, however, are appreciable. For example, in the 
above fit, the color-singlet matrix elements were taken from Buchmiiller-Tye 
wavefunctions at the origin and not considered as free parameters. The over- 
all normalization depends strongly on the charm quark mass, although the 
direct production fraction does not. The short-distance coefficients have been 
expressed in terms of the charm quark mass rather than quarkonium masses. 
This is conceptually preferred by the factorization formalism and numerically 
quite important. 

Even disregarding these uncertainties, the above picture can not be con- 
sidered as complete. Defining a parameter a for the angular distribution as 
before, where 9 is now the angle between the three-momentum vector of the 
positively charged muon and the beam axis in the quarkonium rest frame, one 
finds 0.15 < a < 0.44 for ip' production at y/s = 21.8 GeV and 0.31 < a < 0.63 
for J/ip production at ^/s = 15.3 GeV. [Since the energy dependence is mild, 
these numbers can be used with little error at higher cms energies.] These 
estimates could be made more precise, if the matrix elements (0g^( 3 Po)) and 
(Cg^( 1 S'o)) were individually known. The degree of transverse polarization is 
higher in J/ip production, because the indirect contribution from \ci decays 
yields a purely transversely polarized component to the cross section. These 
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estimates should be compared with the measurement of no visible polarization, 
a s» 0, both for if)' and J/ip. 

At leading order in v 2 , one predicts a % c i : % C 2 production ratio (weighted 
by their branching fractions to decay into J/tp) of 1 : 7, mainly because Xci 
production is suppressed in a s . As for direct J/ijj production, higher order 
corrections in v 2 can be expected to dominate \ci production. Although they 
have not yet been analyzed, they could raise the above ratio to about 1:3, 
still far below the observed ratio (1.4 i-0-4) : 1. Since the Xci ■ Xc2 ratio has 
been measured in only one experimented, an independent confirmation of this 
discrepancy would be welcome. 

The discrepancy between predictions and data on polarization might be 
due to numerically large spin symmetry breaking corrections, which would lead 
to less transverse polarization. One could also suspect that the experiments 
have not yet fully accounted for all systematic errors. In the.E672/E706 exper- 
iment (the only one that provides this piece of informationE3), for instance, the 
acceptance and efficiency varies strongly with angle, while at the same time, 
the acceptance/efficiency curve was determined from a Monte Carlo sample 
of unpolarizcd quarkonia. However, together with the problematic Xci '■ Xc2 
ratio, it seems more likely that the interactions of a color octet cc pair with 
soft gluons as it traverses the target are not understood. Although of higher 
twist in A/m c , such effects must be sizeable as illustrated by the large nuclear 
dependence of the total -0 cross sections, which is dynamically also not un- 
derstood. Since quarkonium formation occurs only after the cc pair left the 
nucleus, it appears as if by the time the quarkonium is formed, the cc pair has 
lost 'memory' that it was produced from two on-shell (transverse) gluons, as 
assumed in leading twist. But even if higher twist effects are sizeable, the po- 
larization problem seems hard to solve, since the coupling of soft gluons to the 
cc pair would be subject to spin symmetry and the counting rules for multipole 
transitions, so that transverse polarization is again hard to avoid. No higher 
twist mechanism has yet been shown to produce predominantly longitudinal 
polarization in the central xf region that dominates the cross section. To clar- 
ify the higher twist nature of the discrepancies, analogous measurements for 
T production would be desirable. As another check, it would be interesting to 
know to what extent Xci is produced with helicity ±2 as predicted at leading 
twist (with small corrections). 

Photoproduction of charmonia in fixed target experiments (or HERA, in 
the low-pt region) is theoretically rather similar to fixed target hadroproduc- 
tion as far as the underlying diagrams and the probed matrix element A$(J/i/)) 
is concerned. A compar.ispn. .of octet mechanisms with data was-jmdertaken 
both in the elastic regionlHJE3'E3 z > 0.95 and the inelastic domaincJ-Ea z < 0.9, 
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where z = p-k^/p-k^ is the energy fraction E^p/E^ in the proton rest frame. In 
the elastic region color-octet contributions are enhanced by ir/a s compared to 
color-singlet contributions, as in fixed target hadroproduction. The extracted 
matrix element As(J/tp) is consistent with the extraction from hadroproduc- 
tion, but substantially smaller than what would have been expected from the 
fit to the Tevatron data in the moderate pt domain discussed earlier. While 
the consistency with hadroproduction is reassuring, the situation with higher 
twist effects is even more delicate in the large- z region than in hadroproduction. 
As octet mechansims, 'diffractive' quarkonium production is also centered at 
z 1 and unsuppressed if not dominant in the endpoint region. In addition, 
the restriction z > 0.95 might not allow the process to be sufficiently inclusive, 
as is necessary in the NRQCD approach. With these remarks of caution in 
mind, the discrepancy in color-octet matrix elements between 'high' energy 
(Tevatron) and 'low' energy fits should probably not be over-interpreted. It is 
likely that energy dependent higher-order radiative corrections (such as small- 
x effects) would, once taken into account, lower the matrix elements extracted 
from the Tevatron fits. Such a reduction would be welcome to. explain a fac- 
tor two discrepancy in comparison with large-pt UA1 data 11-3 and the J/ip 
branching fraction in B decays, see below. 

In the inelastic domain z < 0.9, color-octet processes are parametrically 
suppressed by v 4 , but, due to enhancements in the amplitudes, numerically 
of the same order as color-singlet processes. When z approaches unity, the 
color-octet contributions diverge as expected from a partonic process with a 
lowest order term proportional to 6(1 — z). This divergence indicates that the 
partonic prediction must be considered as a distribution. When smearing is 
applied, together with the constraints on the matrix elements (0g^( 3 Po)) an( i 
(Og^( 1 S'o)) from the elastic peak, the presence of color-octet contributions 
does not appear to be in conflict with the J/ip energy distribution in fixed 
target experiments or at HERA. 

3.3 e + e~ Annihilation and Z° Decay 

Quarkonium production in e + e~ annihilation shows much of the variety of 
production mechanisms in hadro- and photoproduction, while eliminating the 
uncertainties related to hadrons in the initial state. On the other hand, the 
cross sections are not large and there is not much data to compare with. [For 
instance, out of 3.6 million hadronic Z° decays, OPAL has found 24 prompt 
J/1>.] 

The situation when y/s isjipt much larger than 2m c , relevant at CLEO, 
has been investigated in Ref. ell. The parametric dependence of color-octet 
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contributions to the charmonium energy spectrum is analogous to the energy 
spectrum in photoproduction. At leading order in a s , the cc pair can only 
be produced in an octet state, which should be seen as a jet containing a 
charmonium, recoiling against a gluon jet. The energy distribution is peaked 
at the endpoint E max = (s + M^)/(2 v / s), a delta- function smeared over a 
region v 2 E max . Away from the endpoint, color-singlet contributions are more 
important. In addition to the energy spectrum, the J/ip angular distribution 
leads to a striking signature of octet production: Close to the endpoint of the 
energy spectrum, color-singlet production favors J / ip production perpendicular 
to the beam axis, while octet mechanisms produce J/tp predominantly along 
the beam axis. An accurate prediction again requires the values of the matrix 
elements (Oi^( 3 P )) and (O^^Sq)). 

When s 4m 2 , as relevant for Z° decay at LEP, the production patterns 
change. The ratio 4to^/M| is now an important parameter, just as Am 2 . /p 2 at 
largest. The dominant production mechanisms are of fragmentation-type. Al- 
though suppressed by v 4 compared to color-singlet charm fragmentation, gluon 
fragmentation into a color-octet cc pair in a 3 5i state winsEStJ over charm 
fragmentation by a factor of approximately three, depending on the choice of 
values for the color-singlet and octet matrix elements. Gluon fragmentation 
benefits from larger color and flavor factors (all Z° — > qq initiate gluon frag- 
mentation), but most importantly from a double logarithmic enhancement in 
Am 2 /m 2 z from the phase space region, where the fragmenting gluon is collinear 
to and softer than the primary quark. Consequently, color-octet gluon frag- 
mentation could be discriminated from charm fragmentation that yields more 
energetic charmoniaea. 

Present data from LEP E3 have too low statistics to idenitfy a gluon frag- 
mentation contribution unambiguously The central value for the prompt J / %p 
branching fraction Br(Z° -> J/ip prompt + X) = (1.9 ± 0.7 ± 0.5 ± 0.5) ■ 10~ 4 
however, is about a factor of two larger than the sum of all color-singlet contri- 
butions. [Such statements are not quite exact, since the experimental efficen- 
cies and thus results depend on the assumptions on the underlying production 
mechanism. The largest experimental error arises from the determination of 
the prompt production fraction.] The measured energy spectrum is neither 
peaked towards low energies nor towards large energies. Closing the eyes on 
the present errors, one would suspect a roughly equal contribution from charm 
and color-octet gluon fragmentation. 
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3.4 Bottom Decay 

Color-octet mechanisms have been first considered for the inclusive \cj yield 
in B meson decays u. In leading logarithmic approximation, color singlet pro- 
duction is proportional to (2C + — C_) 2 « 0.16, while octet production is 
proportional to (C+ + C_) 2 s» 4.9, where C+ and C_ are related to the Wilson 
coefficients of the current-current operators in the AB — 1 effective Hamilto- 
nian. This numerical enhancement in the Wilson coefficients allows color-octet 
contributions to compete with color-singlet^aatributions, even though the first 
are suppressed by v A in the matrix elementsc3 : Eil. The measured branching frac- 
tion Br(B — > J/ip+X) — (0.80±0.08)% allows us, in principle, to put a bound 
on the least known matrix elements (0%^ ( 3 Pq)) and (Og (^So))- But the 
cancellations in the combination (2C+ — C_) 2 also render its numerical value 
highly sensitive to the choice of scale in the Wilson coefficients. This, together 
with the uncertainties in the wavefunction at the origin, leave an uncertainty 
of about a factor four in the color-singlet contribution. However, even with 
no color-singlet contribution at all and no contribution from (O^^ ( 3 Si)) , the 
Tevatron fittl of (O^^So)) + 3/m 2 (O 8 7/ ^( 3 P )) ~ 6.5 • 1CT 2 would overesti- 
mate the branching fraction by about 50%. Thus, B decays favor the smaller 
values of these matrix elements inferred from fixed target hadro- and photo- 
production. 

Bottomonium decay provides another Jhh .production process. The two 
most important production mechanisms areEirEj T — > \ c j + X, dominated 
by octet \cj production, followed by radiative Xcj decay, and T — > ggg*, 
followed by g* —> cc[ 3 s[ 8 ^}. The second mechanism is similar to color-octet 
gluon fragmentation, although the T mass is not large enough to justify the 
fragmentation approximation. The branching fractions are estimated to be 
7 • 10~ 5 and 2.5 • 10~ 4 , respectively, using the larger value (0%^ ( 3 Si)) — 
0.015 GeV 3 . Their sum is only a factor of two below the upper bound 6.8 ■ 10~ 4 
from ARGUS and consistent with the CLEO result (1.1 ± 0.4) • 10~ 3 within 
2a. 



4 Conclusion 

The magnitude of color-octet vs. color-singlet production cross sections for the 
J/ip production processes discussed in Sect. 3 are summarized in brief in Tab. |[ 
together with the hadronization matrix elements probed in the corresponding 
process. A rather striking observation is that color-octet mechanisms dominate 
every production process. So, why are we discussing them only more than 
twenty years after the discovery of charmonium? 
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Table 2: Importance of color octet contributions in various J/ip production processes and 
the octet matrix elements that could be probed in the corresponding process. 



Process 


Reference 


octet 
singlet 


Matrix elements 


pp, large p t 


Refi 


i- 




30 




pp, moderate p t 


Ref.N 


8 


3 o(8) 3p(8) lo(8) 
1 ' ' J 


Hadroproduction, fixed target 


Ref.yH^ 


2-8 


3p(8) lo(8) 


Photoproduction, z > 0.95 


Ref.Ef 


25 


26 


4 


3p(8) lo(8) 
^0 ' °0 


Photoproduction, z < 0.9 


Ref.^ 


1 


3p(8) lo(8) 


Z° decay 


Ref.Ny 


3 




e + e~ annihilation, large z 


Ref. 


27 




4 


3p(8) lo(8) 


B decay 


Ref.yy 


3 


3c(8) 3p(8) lo(8) 

1 5 ' J 


T decay 


Rcf.yy 


3 


3 5 (8) 



For one thing, many accurate experimental results are rather recent. An- 
other is that the charm mass is not really large enough to make precise predic- 
tions. Thus, a dramatic signature such as in large-p f production at the Teva- 
tron was necessary to induce a reanalysis of other processes. The numbers in 
the table should be considrered as indicative only and are often uncertain by a 
factor of two or more. Even now these uncertainties do not allow us to deter- 
mine all production matrix elements. With the correct theoretical framework 
at hand, further effort, both experimental and theoretical, can be undertaken 
to make the presented quarkonium production picture fully quantitative. 
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